Single spin asymmetries in inclusive hadron production 
from SIDIS to hadronic collisions: universality and phenomenology 
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In a perturbative QCD approach, with inclusion of spin and transverse momentum effects, exper- 
imental data on azimuthal asymmetries observed in polarized semi-inclusive deeply inelastic scat- 
tering and e^e~ annihilations can be used to determine the Sivers, transversity and Collins soft 
functions. By using these functions, within the same scheme, we predict p^p h + X single spin 
asymmetries in remarkable agreement with RHIC experimental data. 

PACS numbers: 13.88.-(-e, 13.60.Hb, 13.85.Ni 



Several azimuthal spin asymmetries have recently been 
measured in different processes: semi-inclusive deeply in- 
elastic scattering (SIDIS) P,[3,[l, inclusive particle pro- 
duction from hadronic collisions [1, IB, 0, B H ITol | 
and from e+e^ annihilations [Tlj . Sizable transverse 
single spin asymmetries (SSAs) and large transverse hy- 
peron polarizations, see e.g. |12|], severely challenged the 
predictions of leading order perturbative QCD (pQCD) 
obtained in coUincar configuration [Tsll, where the con- 
stituent partons and the corresponding observed hadrons 
are assumed to be collinear. Lately, a remarkable success 
in reproducing SSA experimental data has been achieved 
in the framework of non-collinear pQCD, by including 
spin and partonic intrinsic motion effects. This requires 
the introduction of a new class of transverse momen- 
tum dependent (TMD) parton distribution (PDFs) and 
fragmentation functions (FFs). Phcnomenologically, the 
most relevant are, at leading twist, the Sivers [3, 
and Boer-Mulders distribution functions and the 
Collins fragmentation function In fact, in a factor- 
ized pQCD approach, they are responsible for most of the 
azimuthal spin asymmetries observed in SIDIS and Drell- 
Yan (DY) processes, and in e+e" annihilations. Factor- 
ization theorems have been proved for SIDIS and DY 
processes [3], in the kinematical regime where two well 
distinct energy scales are present: a large scale - either 
the photon virtuality or the invariant mass of the lep- 
tonic pair, Q - and a low-moderate scale - the transverse 
momentum, px, of the observed hadron or of the lep- 
tonic pair with respect to the colliding beams - so that 
Aqcd ~ pt < Q- 

Although the large SSAs observed in inclusive sin- 
gle particle production in hadronic collisions motivated 
a huge theoretical effort towards a pQCD treatment of 
these processes with the inclusion of TMD effects, a proof 
of factorization for this case is still lacking. Moreover, it 
has been shown that in a QCD approach the Sivers and 
Boer-Mulders effects vanish unless initial and final state 
interactions between the struck parton and the spec- 
tators are included [H, H^l- These effects, which can 
be accounted for by including appropriate color gauge 



links (Wilson lines) in the invariant definition of the fac- 
torized hadronic correlators, can actually be reabsorbed 
in prefactors associated to the hard partonic cross sec- 
tions. Notice that these prefactors are different for dif- 
ferent processes [2l[. These results, while validating 
the interpretation of the Sivers and Boer-Mulders func- 
tions as the basic mechanisms generating SSAs and az- 
imuthal asymmetries, question their universality prop- 
erties and the whole concept of factorization. In other 
words, for these processes factorization is broken, but 
in a known and predictable way (i.e. through the ac- 
tion of the hard scattering prefactors), so that the to- 
tal polarized and unpolarized cross sections can still be 
calculated in a scheme where factorization should be un- 
derstood at best in a generalized, transverse momentum 
dependent way. Significant progress in this sense has 
been achieved in the context of inclusive almost back- 
to-back two-particle (e.g., dijet, photon-jet) production 
in hadronic collisions [22, 123, [2J|. While these results 
give very useful indications and suggest caution in treat- 
ing single particle production at large pt, they cannot 
give definitive statements on this process, where the low 
energy scale (the intrinsic k±) is integrated over and un- 
observed. 

Alternatively, a more phenomenological approach to 
inclusive single particle production in hadronic collisions 
has been formulated in the context of a generalized par- 
ton model, with the inclusion of spin and transverse mo- 
mentum effects [2^. In this approach factorization 
is assumed as a reasonable starting point, leading twist 
TMD distributions keep their partonic interpretation and 
are therefore expected to be universal and process inde- 
pendent. Although reasonable, these assumptions must 
undergo phenomenological scrutiny by careful compar- 
ison with experimental measurements. In this scheme, 
one should try to reproduce the available data on SSAs 
coming from SIDIS, hadronic collisions and e'^e^ anni- 
hilations with a single set of universal, TMD distribu- 
tions. Historically, this approach was first applied to 



the p^p — > TT 
available at - 



- X process, since low energy results were 
= 20 GeV from the E704 CoUabora- 
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tion. Later on, high energy regimes became available 
at RHIC y/s = 200 GeV experiments. Recently, the 
same approach has been independently applied to the 
Sivers and Collins azimuthal asymmetries measured by 
the HERMES and COMPASS Collaborations, and to the 
double hadron azimuthal correlations measured in e^e~ 
annihilations. High quality global fits to all experimen- 
tal data have allowed the simultaneous determination of 
the Sivers [l^l and transversity distributions and of the 
Collins fragmentation function [28j . 

In this paper, for the first time, we investigate the 
possibility of reproducing unpolarized cross sections and 



SSAs in proton-proton collisions, by using the Sivers and 
transversity distribution functions, and the Collins frag- 
mentation function, as determined by fitting SIDIS and 
e+e~ experimental data, without any modification. No 
fit to pp collision data will be used nor performed, so that 
all our results may be considered as genuine predictions. 

Let us briefly sketch the basic relations adopted in our 
approach. For a detailed and complete treatment we refer 
to the original papers [1^, [2^ . The invariant differential 
cross section for the polarized process p^'^p ^ tt + X can 
be written as 
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The SSA can be expressed as the ratio = {da^ — da^)/{da^ + da^). In Ref. [l^ it was shown that in the 
TMD generalized parton model several terms, other than those present in the coUinear configuration, contribute to 
both the numerator and the denominator of Ajf. These contributions were neglected in earlier TMD parton models 
which lacked proper inclusion of spin effects. However, Refs. [l^, [29l| have shown that for both the E704 and RHIC 
kinematics, the numerator of the SSA is dominated by the Sivers effect, with a small or negligible contribution from 
the Collins effect, while the denominator is essentially given by the TMD unpolarized contribution which generalizes 
the usual coUinear one. Therefore, in the case of interest here, Eq. ([T]) takes the simple form 
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Notice that the contribution due to the convolution of the 
transversity distribution with the Collins fragmentation 
function [26|, has actually been included in our calcula- 
tions. However, this effect is marginal in our results, with 
the exception of the BRAHMS An of Fig. H AU other 
terms are totally negligible in the kinematical configura- 
tions considered, as their azimuthal phase factors become 
very complex and in general, under full k± integrations, 
wash out the corresponding contributions. Consequently, 
the dominant terms at numerator and denominator are 
the Sivers effect and the unpolarized term, which have 
the simplest azimuthal phases. 

In order to perform the calculations, we strictly fol- 
low Refs. [23, [ill where the azimuthal asymmetries mea- 
sured by the HERMES and COMPASS Collaborations 
for semi-inclusive pion and charged hadron production 
have been used to extract the Sivers function and, via 
a combined analysis of Belle data on e~^e~ —> tttt + X, 
the transversity distribution and the Collins fragmenta- 
tion function. It is important to recall that in SIDIS the 
Sivers and Collins contributions to the asymmetry can 



be disentangled by considering suitable moments of the 
azimuthal distributions of the observed hadrons. 

In what follows, we will adopt the same Sivers and 
transversity distributions and the Collins FF as in 
Refs. 27, 28]. Our aim is to calculate the SSA for pion 
and kaon production in polarized pp collisions, for the 
RHIC kinematical range. Notice that no fit is performed, 
we will simply compare our predictions to data on un- 
polarized cross sections and SSAs from the STAR and 
BRAHMS Collaborations. AU TMD functions have a 
simple ansatz form in which the k± and lightcone fraction 
dependences are factorized, with a flavor independent 
Gaussian shape for the k±^ term. For unpolarized PDFs 
and FFs the Gaussian widths are flxed to (fc^) = 0.25 
(GeV/c)2 and {pD = 0.2 (GeV/c)^, respectively [H . 
As in Ref. [13, we adopt the MRSTOl PDF set [3l| and 
the Kretzer FF set [13]. Following Ref. [2^ we fix the fac- 
torization scale to /i = Pct/2, where PcT is the transverse 
momentum of the fragmenting parton c in the partonic 
cm. frame. Notice that, although our calculations are 
performed at lowest order in pQCD, no additional K- 



3 



factors are introduced. 

In Fig. [1] panel (a), we show the unpolarized cross 
section for neutral pion production in pp collisions, at 
^/s = 200 GeV and two different average pseudorapidi- 
ties, (77) — 3.3 and {rj) — 3.8, as a function of E.^. The 
corresponding SSAs, A]y{p^p tt^ + X), are plotted in 
panels (b), (c) and (d) as a function of xp at fixed ra- 
pidity values, (b), and as a function of pr at different 
bins in xp, (c) and (d). The curves are our predictions, 
compared to STAR experimental data [l^, [s^ • In this 
case the Collins effect is totally negligible. 

In Fig. [21 (a) and (b) panels, we compare our pre- 
dictions to BRAHMS data on unpolarized cross sections 
for charged pion production [sj , as a function of pr at 
y/s = 200 GeV and two different rapidities, y — 2.95 and 
y — 3.3. On the (c) and (d) panels, we plot the anal- 
ogous results and data for An [9], as a function oi xp 
and for two different cm. scattering angles, 6 — 2.3° and 
6 — 4.0°. In particular, the total contribution obtained 
by adding the Sivers and the Collins effect (thick lines) is 
shown together with the Collins contribution alone (thin 
lines). 

Finally, in Fig. [3] we compare our predictions for un- 
polarized cross sections, (a) and (b), and SSAs, (c), to 
the BRAHMS data for charged kaon production [3, l3J] . 
Notice that for Ajv we consider only the Sivers effect as 
the Collins contribution is small and the unknown kaon 
Collins functions are still under study. 

Let us comment on these results: 

1) On the whole, with the exception of low px STAR 
data, our predictions are in remarkably good agreement, 
both in size and sign, with RHIC data. In particular, 
there are no evident contradictions that could be inter- 
preted as signals of sizable factorization and universality 
breaking between SIDIS and pp collisions. We believe 
this information may be very useful in the context of 
the recent theoretical developments: in fact, although 
factorization-breaking terms could in principle be there, 
it is not easy to infer how sizable and significant their 
contributions should be for the processes and kinemati- 
cal configurations where data are presently available. 

2) Although remarkable, one should not overestimate the 
significance of our results. First of all, the HERMES and 
COMPASS data on SIDIS used in the fits of Refs. [H,!!! 
are collected mainly in the region of low- intermediate val- 
ues of the Bjorken variable {xb < 0.3). Therefore, the 
fits strictly constrain our parameterization for the Sivers 
function only in the low x region. Consequently, with 
these parameterizations we are unable to reproduce the 
low energy results on the SSA for charged pions at large 
Xp oi the E704 collaboration. Similarly, with the average 
k±'s for PDFs and FFs extracted from SIDIS, we under- 
estimate the corresponding unpolarized cross sections by 
a factor which cannot simply be attributed to pQCD K- 
factors. Further advances in this direction would require 
a simultaneous fit of SIDIS and pp data, which would 
give tighter constraints on the large x behaviour of our 



parameterizations. 

3) In the color gauge invariant TMD approach of Ref. [2lj 
gauge links result in numerical prefactors which can be 
absorbed in generalized partonic cross sections. These 
factors may be different for different processes, therefore 
breaking factorization. No calculations have been per- 
formed within this scheme for the single inclusive parti- 
cle production in hadronic collisions, but only for the 
double inclusive case. However, it is possible that in 
the kinematical configurations considered here, the dom- 
inant partonic contributions are those with positive pref- 
actors, close to one. Therefore, our findings cannot be 
interpreted as contradicting the color gauge invariant ap- 
proach; rather, they could indicate that in these situa- 
tions the generalized parton model and the color gauge 
invariant approach give very similar results. Unambigu- 
ously different predictions from the two approaches could 
instead be found in SSAs in DY process and in photon- 
jet production in pp collisions, see e.g. Refs. [27, 35, 36J. 

4) For the pp case (and for SIDIS at large transverse 
momentum of the observed hadron, comparable with Q) 
there is a well known alternative approach, the coUinear 
twist-three formalism [13]. Ref. [Ig showed that results 
comparable or better than those found in a TMD ap- 
proach can be obtained for SSAs in pp collisions. Let us 
point out, however, that in Ref. [s^ all available pp data 
were included in the fit, while our results can be consid- 
ered as genuine predictions. Concerning the low energy 
E704 data, this approach faces similar problems both for 
SSAs and unpolarized cross sections. 

5) The fit of SIDIS data used here was performed before 
data on Sivers and Collins asymmetries for neutral pions 
and kaons were made available by the HERMES collab- 
oration. Also COMPASS results on separated pions and 
kaons became available afterwards. This could partially 
explain the poorer agreement with kaon production data 
from BRAHMS. Nevertheless, let us stress that kaon data 
are much more challenging in SIDIS itself, since they are 
crucially dependent on the strange and m-, d-sea quark 
contributions to the Sivers function. Notice that non- 
leading contributions in the fragmentation functions are 
presently known with less accuracy also in the coUinear 
unpolarized case. Some progress in this direction has 
been recently achieved and a new set of fragmentation 
functions is available for kaons [39| which differ substan- 
tially from those adopted in this paper. A detailed study 
of these effects is beyond the aim of this paper and re- 
quires a new fitting procedure and an updated analysis, 
which will be discussed elsewhere [i^l . 

In conclusion, this phenomenological analysis shows 
that most of the available high energy data on unpolar- 
ized cross sections and SSAs for single inclusive particle 
production in SIDIS and in pp collisions can be well re- 
produced in the framework of a generalized parton model, 
with inclusion of spin effects and leading twist TMD dis- 
tributions, in particular the Sivers function. 
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Figure 1: Panel (a): Invariant differential cross section for pp ^ tv'^ + X at y/s — 200 GeV and two pseudorapidity values, 
rj = 3.3 and r) = 3.8, vs. Et^- Data are from STAR [s^]. Curves are obtained adopting the unpolarized fcx-dependent PDFs 
and FFs of Ref. !30|]. Panels (b), (c) and (d): An for pp — > vr'^ + X at y/s — 200 GeV and two pseudorapidity values, 77 = 3.3 
and 77 = 3.7, vs. xf (b), and at different xp values vs. pr, (c) and (d). Data are from STAR [lO|]. Curves are obtained using 
the Sivers functions as determined in Ref. [23| by fitting SIDIS data. 
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Figure 2: Panels (a) and (b): Invariant differential cross section for pp tt* + X at = 200 GeV and two pseudorapidity 
values, y = 2.95 and y = 3.3, vs. pr. Data are from BRAHMS M. Curves are obtained adopting the same choices as in 
Fig. 1(a). Panels (c) and (d): Am for pp tt* +X at = 200 GeV and two different scattering angles, 6 = 2.3° and 6 = 4°, 
vs. Xf- Data are from BRAHMS .9]. Thick curves are obtained adding the Sivers effect, as extracted from SIDIS data in 
Ref. [27|. and the Collins effect coupled with the transversity function, as extracted from a global fit of SIDIS and e^e" data 
in Ref. [28| |. The thin lines show the Collins effect: notice that its sign is opposite w.r.t. the Sivers contribution. 




Figure 3: Panels (a) and (b): Invariant differential cross section for pp + X ai ^/s = 200 GeV and two pseudorapidity 

values, y = 2.95 and y = 3.3, vs. pr- Data are from BRAHMS [3J|. Curves are obtained adopting the same choices as in 
Fig. 1(a). Panel (c): Am for pp — > + X at y/s = 200 GeV and fixed scattering angle, 9 = 2.3°, vs. xf- Data are from 
BRAHMS i9|. Curves are obtained with the Sivers function as determined in Ref. |27j l by fitting SIDIS data. 
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We remark that while SIDIS data on pion and charged 
hadron production have been used for the fits, it is not so 
for the pp case, where all our results are genuine predic- 
tions. Their agreement with data seems to indicate that 
the role of possible factorization-breaking terms may be 
marginal for the processes and kinematical ranges con- 
sidered here. We believe that this phenomenological in- 
formation may presently be very useful, given the rapid 
development of this field of research. Our results might 
indicate that our approach and the pQCD color gauge- 
link formalism cannot be distinguished on the basis of 
these processes. One should rather look at Sivers sin- 
gle spin asymmetries in DY processes, in prompt photon 
or, as suggested in Ref. |36|, photon-jet production in pp 
collisions. These reactions might also be useful in dis- 
entangling our approach from the coUinear twist-three 
formalism. 

In future, data at large xf{xb) in collisions(SIDIS) 



would be very helpful in further constraining our param- 
eterizations, testing our model more severely. Moreover, 
new kaon data challenge the phenomenology of SSAs, as 
they definitely involve the role of sea quark TMD distri- 
butions. This issue will be studied in Ref. [i^]. 
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